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Abstract. Monte Carlo simulation is an essential tool in emission tomography 
that can assist in the design of new medical imaging devices, the optimization of 
acquisition protocols, and the development or assessment of image reconstruction 
algorithms and correction techniques. GATE, the GeantJf Application for 
Tomographic Emission, encapsulates the Geant4 libraries to achieve a modular, 
versatile, scripted simulation toolkit adapted to the field of nuclear medicine. 
In particular, GATE allows the description of time-dependent phenomena such 
as source or detector movement, and source decay kinetics. This feature makes 
it possible to simulate time curves under realistic acquisition conditions and to 
test dynamic reconstruction algorithms. This paper gives a detailed description 
of the design and development of GATE by the OpenGATE collaboration, 
whose continuing objective is to improve, document, and validate GATE by 
simulating commercially available imaging systems for PET and SPECT. Large 
effort is also invested in the ability and the flexibility to model novel detection 
systems or systems still under design. A public release of GATE licensed under 
the GNU Lesser General Public License can be downloaded at the address 
lirttp:// www-lphe.epfi.ch/GATE/ Two benchmarks developed for PET and 
SPECT to test the installation of GATE and to serve as a tutorial for the 
users are presented. Extensive validation of the GATE simulation platform has 
been started, comparing simulations and measurements on commercially available 
acquisition systems. References to those results are listed. The future prospects 
toward the gridification of GATE and its extension to other domains such as 
dosimetry are also discussed. 
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1. Introduction 

Emission tomography is becoming increasingly important in modern medicine for 
both diagnostic and treatment monitoring, with a demand for higher imaging quahty, 
accuracy, and speed. Recently enhanced by the wider availability of powerful computer 
clusters, Monte Carlo simulations have become an essential tool for current and future 
emission tomography development. Examples of research areas benefiting from these 
developments are the design of new medical imaging devices, the optimization of 
acquisition protocols, and the development and assessment of image reconstruction 
algorithms and correction techniques. 

Currently there are numerous Monte Carlo simulation packages for either PET 
(Positron Emission Tomography) or SPECT (Single Photon Emission Computer 
Tomography), each with different advantages, disadvantages, and levels of reliability 
(Buvat and Castiglioni 2002). Accurate and versatile general-purpose simulation 
packages such as GeantS (Brun et al 1987), EGS4 (Bielajew et al 1994), MCNP 
(Briesmeister 1993), and most recently Geant4 (Agnostelli et al 2003) are available. 
These packages include well-validated physics models, geometry modeling tools, and 
efhcient visualization utilities. However, it is quite difficult to tailor these packages 
to PET and SPECT. On the other hand, the dedicated Monte Carlo codes developed 
for PET and SPECT suffer from a variety of drawbacks and limitations in terms of 
validation, accuracy and support (Buvat and Castiglioni 2002). While an adaptation 
of EGS4 to radiation therapy applications exists (Rogers et al 1995, Kawrakov and 
Rogers 2003), there are no dedicated PET or SPECT Monte Carlo programs that are 
detailed and flexible enough for realistic simulations of emission tomography detector 
geometries. SimSET (Harrison et al 1993), one of the most powerful dedicated codes 
enabling PET and SPECT simulations, enables a precise and efficient modelling 
of physics phenomena and basic detector designs (e.g. ring detectors and planar 
detectors), but it has limitations with respect to the range of detector geometries that 
can be modeled. For example, a detector ring can not bo subdivided into individual 
crystals and the gaps between the crystals and the grouping of crystals into blocks 
cannot be modeled. In addition, neither SimSET, nor any other publicly available 
codes accoLmt for time explicitly, which limits their use for modeling time dependent 
processes such as tracer kinetics or bed motion. 

Clearly, a Monte Carlo code capable of accomodating complex scanner geometry 
and imaging configurations in a user-friendly way, while retaining the comprehensive 
physics modeling abilities of the general purpose codes is needed. Furthermore, the 
need is to have a platform that can model decay kinetics, deadtime, and movement, 
while benefiting from the same versatility and support as that of the general-purpose 
simulation codes. Object-oriented technology appeared to be the best choice to ensure 
high modularity and re-usability for a PET and SPECT simulation tool. Therefore, 
we selected the simulation toolkit developed in C+-f by the Geant4 collaboration 
(Agnostelli et al 2003), and decided to foster long-term support and maintenance 
by sharing code development among many research groups forming the OpenGATE 
collaboration. 

This paper presents a detailed description of the design and development of a 
Monte Carlo tool by the OpenGATE collaboration which satisfies the requirements 
mentioned above. It was launched at first (Strul 2001, Strul 2001b) as an aid in 
the design of the ClearPET prototype scanners being developed by the Crystal Clear 
collaboration (Ziemons et al 2003). GATE, the Geant4 Application for Tomographic 
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Emission (Santin et al 2003, Strul et al 2003, Assie et al 2004), incorporates the 
Gcant4 libraries in a modular, versatile, and scripted simulation toolkit that is 
specifically adapted to the field of nuclear medicine. A public release of GATE licensed 
under the GNU Lesser General Public License {LGPL n.d.) can be downloaded at the 
address http://www-lphe.epfl.ch/GATE/. 

2. GATE basics 

GATE was designed with several objectives in mind. First, the use of the 
GATE software should not require any knowledge of C+-I-. End-users from the 
nuclear medicine community should be able to use GATE without worrying about 
the programming details. Second, as many nuclear medicine diagnostic techniques 
share similar concepts, GATE software components should be general enough to be 
reused from one context to another. Last, GATE should be modular, and thus be able 
to evolve as new applications are envisioned. 

2.1. Architecture 

The requirements discussed above are met using a layered architecture. The core 
of GATE, developed in C++, defines the main tools and features of GATE. The 
application layer is an extensible set of CH — h classes based on the GATE core. On top 
of the application layer is the user layer, where end-users can simulate experiments 
using an extended version of the Geant4 scripting language. 




Developer 
layer 



Figure 1 . Sketch of the layered architecture of GATE. 

The GATE developer layer consists of the core layer and the application layer. It 
is built from the various classes that provide the most general features of GATE. These 
classes define which tools are available, what developers can do, and how they can do 
it. The core layer includes some base classes that are common or even mandatory in all 
Geant4-based simulations, such as those involved in the construction of the geometry, 
the interaction physics, the event generation, and the visualization management. In 
addition, the core layer includes classes that are specific to GATE simulations, such as 
the GATE virtual clock for time management. Thus, the core layer defines the basic 
mechanisms available in GATE for geometry definition, time management, source 
definition, detector electronics modelling, and data output. 

The application layer is composed of classes derived from the base classes of the 
core layer to model specific objects or processes. For example, the core layer defines 
the base class for volumes, and the application layer comprises all the derived classes 
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for modeling specific volumes, including boxes, spheres, cylinders, and trapezoids. 
Similarly, the application layer includes all the specific movement models derived from 
the movement base class, including translations, rotations, orbits, and oscillations. 
Thus, the range of features available in GATE can increase as new application classes 
are developed, while the general structure remains unaS'ected. 

In the user layer, Geant4 provides mechanisms for running simulations both 
interactively or batch-wise using scripts. An important principle of GATE is that 
each class must provide dedicated extensions to the command interpreter class, so 
that the functionality provided by the class is available through script commands. 
The end-users of GATE therefore do not have to perform any C++ coding. The 
complete set-up of a nuclear medicine experiment can be easily defined using the 
script language, as shown in Figure El 




# CRYSTAL 

/t lock/daugters/name crystal 
/b lock/daugters/insert box 
/crystal/geometry /setXLength 30 mm 

/crystal/geometry /set YLength 4.4 mm 

/crystal/geometry/setZLength 4.4 mm 

/crystal/setMateiial BGO 

# REPEAT CRYSTAL INSIDE BLOCK 
/c rystal/repeaters/insert cubic AiTay 
/crystal/cubicArray/setRepeatNumberX 1 

/crystal/cubicArray/setRepeatNumbei'Y 8 

/c rystal/cubicAnay/setRepeatNumbeiZ 8 

/crystal/cubicAnay/setRepeatVector 0. 4.5 4.5 mm 

# REPEAT BLOCK INSIDE RING 
/block/placement/setTranslation 125. 0. 0. mm 
/block/repeaters/insert ring 
Zblock/ring/setRepeatNumber 1 8 



Figure 2. Ring detector geometry obtained by using the command scripts 
displayed in the figure. This script models a detector-block, where 4.4 X 4.4 mm 
BGO crystals of 30 mm thickness are repeated in an 8 X 8 matrix. A cubic array 
repeater command produces an 8 X 8 cubic array of crystals from a single crystal 
with a crystal picth of 4.5 mm in both directions. A ring repeater places the block 
at 125.0 mm from the origin on the X-axis and replicates it eighteen times about 
a cylinder around the Z-axis. 



2.2. Systems 

When defining the geometry for a tomograph, specific guidelines with respect to 
the geometrical hierarchy of the tomograph components must be followed, so that 
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the Geant4 particles interaction histories, called hits, occuring in the detector, can 
be processed to realistically mimic detector output. Most PET scanners are built 
following comparable concepts: one or more rings, each ring consisting of several 
scintillator blocks, each block being subdivided in crystal pixels. For SPECT, similar 
concepts exist: a gamma camera with a continuous or pixelated crystal, and a 
collimator. Most of these geometrical concepts are common to many different imaging 
systems. To facilitate the hierarchical description of a tomograph, predefined global 
systems are used. A system is defined as a family of geometries compatible with 
one or several data output formats. The main property of a system is that its 
geometry description is supported by specific list-mode or histogrammed (sinogram 
or projection) data output formats. Currently there are five systems available in 
GATE: one for SPECT, three for PET — two for block detector geometries and one 
for continuous pixellated geometries — and a generic system appropriate to model 
novel tomographic paradigms. This latter system is completely open and provides 
only basic building blocks for the definition of a tomographic experiment. 

2.3. Management of time and movements 

One of the distinctive features of GATE is the management of time-dependent 
phenomena (Santin et al 2003, Strul et al 2003). The synchronization of the source 
kinetics with the movement of the geometry thus allows for the simulation of realistic 
acquisition conditions including patient movement, respiratory and cardiac motions, 
scanner rotation, or changes in activity distribution over time. Dealing with time 
in GATE includes: (a) defining the movements associated with the physical volumes 
that describe the detector and phantom; (b) describing the radioactive sources; and 
(c) specifying the start and stop times of the acquisition (which are equivalent to the 
start and stop times in a real experiment). 

The Geant4 geometry architecture requires the geometry to be static during 
a simulation. However, the typical duration of a single event is very short when 
compared to movements in the geometry model or bio-kinetics. Movements are 
synchronized with the evolution of the source activities by subdividing the acquisition 
time frames (typically of the order of minutes or hours) into smaller time steps. At 
the beginning of each time step, the geometry is updated according to the requested 
movements. During each time step, the geometry is held at rest and the simulation 
of the particle transport proceeds. Within the time steps, the source is allowed to 
decay so that the number of events decreases exponentially from one time step to 
the next, and within the time steps themselves. The proper timing of the simulated 
event sequence is a key feature for modeling time-dependent processes such as count 
rates, random coincidences, event pile-up, and detector deadtime (Simon et al 2004). 
Between time steps, the position and the orientation of a subset of daughter volumes 
can be changed to mimic a movement such as a rotation or a translation. These 
displacements are parametrized by their velocity. It is the responsibility of the user to 
set the time step duration short enough to produce smooth changes. Combinations of 
translations and rotations allow the simulation of complex acquisition trajectories of 
the detectors such as parameterized eccentric rotations. 

Geant4 does not allow the movement of sources. Therefore, in GATE, an emission 
volume is defined so that it encompasses the actual source's volume throughout its 
range of displacement. To enable movement of the activity distribution, an additional 
volume is defined to confine the emission within the intersection of the emission and 



GATE: a simulation toolkit for PET and SPECT 



7 



the confinement volumes. This confinement volume defines the shape and size of the 
actual source and moves within the emission volume. 

3. Physics 

3.1. Radioactive sources 

A source in GATE is defined by its particle type {e.g. radionuclide, gamma, positrons, 
etc.), position (volume), direction (solid angle), energy (spectrum), and activity. The 
lifetime of a radioactive source is usually obtained from the Geant4 database, but it 
can also be set by the user to approximate a decay source through the emission of its 
decay products {e.g. positrons or gammas). 

The activity determines the decay rate for a given source during the simulated 
acquisition time. Radioactive decay of radionuclides with secondary particle emission 
is performed by the Geant4 Radioactive Decay Module (RDM), which has been 
modified so that GATE source manager maintains control over the definition of decay 
time. Continuous event time flow is obtained by using a virtual clock that defines 
an absolute time t used to initialize the Geant4 internal tracking time. Random 
time intervals 5t between events that occur at time t are sampled from a exponential 
distribution: 

p{5t) = A{t) exp(-yl(i) 5t) (1) 

where A{t) = Aoexp((t — to) / r) is the source activity at time t, Aq is the user 
defined source activity at time to and r the lifetime. When the resulting decay time 
exceeds the end of the current time step, the run is terminated and a new one is started, 
allowing for the synchronization of the sources with the geometry movements. Multiple 
sources can be defined with independent properties. For each event, a proposed time 
interval is sampled for each source according to Eq. ^ and the shortest one is chosen 
for the primary generation. The overall behavior of this mechanism is such that, 
at all times, the relative importance of each source is proportional to its activity, 
while the overall time interval sampling is determined by the total activity of all 
sources. Voxelized phantom or patient data can be used as sources to reproduce 
realistic acquisitions: emission data are converted into activity levels, and GATE can 
read in voxelized attenuation map and converts the gray scale into material definitions 
using an analogous translator. 

3.2. Positron emission 

GATE includes 2 modules dedicated to PET (Jan 2002). The first uses the von 
Neumann algorithm (von Neumann 1951) to randomly generate the positron energy 
according to the measured spectra. This method greatly increases the speed of the 
simulation by bypassing the decay of radionuclides process used by Geant4. The /3+ 
spectra of 3 commonly used radionuclides (^^C,^^0, and ^^F) have been parametrized 
in GATE according to the Landolt-Bornstein tables (Behrens and Janecke 1969). 

The second module deals with the acollinearity of the two annihilation photons, 
which is not accounted for in Geant4 . In GATE, acollinearity is modelled using a 
0.58° full width at half maximum (FWHM) Gaussian blur. This width corresponds 
to experimental values measured in water (Iwata et al 1997). 
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3. 3. Interaction modelling with standard energy and low energy packages 

A material database file contains all material parameters required by Geant4 to 
calculate the interaction cross-sections and is easily modified by the user. In contrast 
to Geant4, GATE only uses natural isotopic abundances. The fact that these cannot be 
modified by the user has little bearing on GATE applications since isotopic abundances 
are unimportant in low to mid-energy photon and charged particle interactions. 
The electromagnetic interactions used in GATE are derived from Geant4 . The 
electromagnetic physics package manages electrons, positrons, 7-rays, X-rays, optical 
photons, muons, hadrons, and nuclei. As in Geant4, GATE can use two different 
packages to simulate electromagnetic processes: the standard energy package, and the 
low energy package. In the standard energy package, photoelectric effect and Compton 
scatter can be simulated at energies above 10 keV. Under 100 keV however, relative 
errors on the cross-sections are higher than 10% and can raise up above 50% (Jan 2002 
and Lazaro 2003). The low energy package models photon and electron interactions 
down to 250 eV and includes Rayleigh scattering. For biomedical applications, it 
provides more accurate models for the electromagnetic interactions. However, this 
comes at the price of increased computing time. 

3.4- Secondary production cuts 

GATE inherits the Geant4 capability to set thresholds for the production of secondary 
electrons, X-rays and delta-rays (Agnostelli et al 2003). In biomedical applications, 
eliminating the secondary particles whose initial energy is below the production 
threshold increases the computing efficiency. 

Because low energy processes generate more secondary particles than standard 
energy processes, cuts affect simulation speed more strongly when applied with the 
low energy package. Turning off the production of electrons. X-rays, and delta-rays 
by setting high thresholds may result in a substantial increase in computing speed for 
a typical simulation of a PET scanner. In many cases, the accuracy of the simulation 
at the level of single or coincidence photon counting is preserved. 

4. Digitization 

Digitization is the process of simulating the electronics response of a detector within a 
scanner. This involves the conversion of the charged particle and photon interactions 
into energy bins, detection positions, and coincidences. To do this, portions of the 
scanner geometry are designated as sensitive detectors, which record interactions 
within these regions. The digitizer chain then processes these recorded interactions 
and produces counts and coincidences. The sensitive detectors and digitizer chain are 
described below. 

Sensitive detectors 

Sensitive detectors are used to store information about particle interactions (hereafter 
referred to as hits) within volumes. GATE only stores hits for those volumes that have 
a sensitive detector attached. Two types of sensitive detectors are defined in GATE: 
the crystal sensitive detector {crystalSD) is used to generate hits from interactions 
that occur inside the detector portions of the scanner (Figure Ol. The phantom 




Figure 3. Particle interactions in a crystalSD attached to a scintillator block and 
a phantomSD attached to a volume filled with tissue material. The trajectory of 
particle A shows 1 hit in the phantomSD (Hit al) and 4 hits in the crystalSD 
(Hit a2 to Hit a5). Particle B does not interact within a sensitive volume, thus 
no hit information is stored. 



sensitive detector {phantomSD) is used to detect and count the Compton and Rayleigh 
interactions occurring within the scanner's field-of-view (FOV). 

4-2. Digitizer chain 

The digitizer chain mimics a reahstic detection process by building the physical 
observables from the hits. The observables of each detected event are energy, position, 
and time of the interaction. The digitizer consists of a chain of processing modules 
that takes a list of hits from the sensitive detectors and transforms them into pulses 
referred to as singles. The key elements of this chain are now briefly described. 

4-2.1. Hit adder A particle entering a sensitive detector can generate multiple hits, 
as shown in Figure 13 For instance, a gamma ray interacting within a scintillation 
crystal can generate hits corresponding to several Compton scattering events and a 
photoelectric absorption. The hit adder sums the deposited energy of these hits within 
a sensitive detector to yield a pulse. The position of the pulse is calculated from the 
energy- weighted centroid of the hit positions, and the time of the pulse is set to that 
of the first hit within the volume. If a particle interacts in several sensitive detectors 
the hit adder will generate a list of pulses, one for each sensitive detector. 

4.2.2. Pulse reader With the exception of one-to-one readout schemes, scanners often 
have a readout segmentation different from the detector segmentation. To simulate 
this, a pulse reader module adds the pulses together within a user-defined group of 
sensitive detectors. This yields a pulse containing the total energy deposited within 
the group of sensitive detectors. The position of this pulse is set to that of the pulse 
from the adder that has the largest energy (winner-takes-all paradigm) . 



GATE: a simulation toolkit for PET and SPECT 



10 



4-2.3. User modules Following the hit adder and the pulse reader, which regroup 

the hits into pulses and then sum pulses, the remaining modules of the digitizer chain 
transform these pulses into the physical observables of the scanner {i.e. singles). These 
modules are discussed below. 

Energy resolution: This module applies a Gaussian blur to the energy of the pulse, 
E, with an FWHM of Rx E. The FWHM energy resolution, R, is given by: 

R=?^ (2) 



where Rq is the user defined FWHM energy resolution at a given energy, Eq. 

A more elaborate model propagates the relative variances of the physical processes 
involved with light collection and detection in a spectrometric chain as: 

^^ = "^+iV + iV;^ + iV^ (3) 
where ve^ vj^, w,,, v^, and vm are the relative variances on E, on the number 
of scintillation photons N, on the light collection efficiency rj, on the quantum 
efficiency of the photo-detector e, and on the gain of the photo-detector M. Both 
the light collection efficiency and the quantum efficiency are binomial processes 
with probabilities of rj and e, respectively. If the scintillation process is assumed 
to follow a Poisson law with a mean equal to N — E x LY, where LY is the light 
yield of the scintillator, then relative variance on E is given by: 

In the case of a photomultiplier, vm is approximately 0.1 and the FWHM energy 
resolution is estimated using the equation: 



where Ri is the intrinsic resolution of the scintillator (Kimtncr et al 2002). 

Energy window: Upper and lower energy thresholds can be set for several energy 
windows by using multiple processor chains. These thresholds are applied using 
either a step or sigmoid function. 

Spatial resolution: For SPECT, spatial resolution is modeled using a Gaussian blur 
of the position. For PET, interaction position is calculated by the pulse reader 
which simulates the intrinsic spatial resolution of the detector. More elaborate 
models can be derived for continuous crystal detector PET systems (Staelens et 
al 2004). These models are currently under development. 

Time resolution: Simulation of time jitter can be obtained using a Gaussian blur 
of the pulse time. 

Deadtime: Both paralyzable and non-paralyzable deadtimes can be modeled 
explicitely on an event-by-event basis. While these models represent the idealized 
behavior, they correctly predict the theoretical lifetimes for both types of 

deadtimes (Simon et al 2004). 

Coincidence: At the end of a digitizer chain a coincidence sort can be added to 
find pairs of singles that are in coincidence. Pairs of singles can be considered 
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coincidences whenever the time interval between the singles is less than a user- 
defined coincidence window. Each single is stored with its corresponding event 
number. If the event numbers of the singles associated in a coincidence are 
different, this is a random coincidence. A similar flag exists for Compton 
scattered events. The Compton scatter flag can be used to differentiate true from 
scattered coincidence pairs that have identical event flags. Multiple coincidences 
corresponding to more than 2 singles within the same coincidence window are 
discarded. 

4.3. DigiGATE 

In GATE standard operation mode, primary particles are generated by a source 
manager and then propagated through the attenuating geometry before generating 
hits in the sensitive detectors, which are then processed by the digitizer chain. While 
this operating mode is suitable for many purposes, it is inefficient for the optimization 
of the digitizer chain parameters. This is best done by comparing the results from 
different sets of digitizer parameters using the same series of hits. To perform this 
specific task, GATE offers an operating mode named DigiGATE. In this mode, hits 
are read from a data file generated by GATE and fed directly into the digitizer 
chain. The same command scripts are used for both the hit generation simulation 
and DigiGATE simulations. Thus, all conditions are kept identical in the simulations 
including time-dependence. 

5. Simulation benclimarks 

Two benchmarks, one for PET and one for SPECT, are included in the 
GATE distribution. These benchmarks check the integrity of the installation or 
upgrade, and also allow for the comparison of CPU performance on different computing 
platforms. In addition, they provide examples of how to use the main features of 
GATE to simulate PET or SPECT experiments. Furthermore, they serve as examples 
on how to analyze output. Each benchmark consists of macros to run the simulation, 
analyze simulation output, and generate figures. In addition, a set of baseline figures 
are included for a comparison between the user's results with those from a correct run. 

5.1. PET benchmark 

The PET benchmark (Figure 0} simulates a whole-body scanner that does not 
correspond to any existing system. Rather, it serves as a simple system that contains 
the majority of GATE features. It consists of eight detector heads arranged within 
a 88 cm diameter by 40 cm axial length octagonal cylinder. Each head is made 
of 400 detector blocks and each block is a 5 x 5 array of dual-layer LSO-BGO 
crystals. The heads are equipped with partial septa that rotate in a step-and-shoot 
mode. The phantom in this benchmark is a 70 cm long water cylinder with one ^^F 
(half-life = 109.8 min) and one ^^O (half-life — 2.03 min) line source each with an 
activity of 100 kBq. The simulated acquisition is 4 min in duration, which represents 
approximately two ^^O half- lives. The source activities are set such that the PET 
benchmark will run in about 12 CPU hours on a 1 GHz processor. The acquisition is 
divided into two 2 min frames; after the first frame, the gantry rotates by 22.5 deg. 
Only coincident events are recorded, using a coincidence time window of 120 ns. This 
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Figure 4. Illustration of the PET benchmark setup. 



large window is used in order to record a large number of random coincidences. The 
lower and upper energy thresholds are set to 350 and 650 keV, respectively. 

The standard electromagnetic package of Geant4 is used in this benchmark. Only 
the Rayleigh interactions are modeled using the low energy package. To speed up the 
simulation, X-rays and secondary electrons are not tracked. 

Approximately 3.7 x 10'' decays occur during the simulated acquisition and around 
7.0 X 10^ coincidences are recorded and written in a ROOTfAc (Brun and Radcmakers 
1997). Based on the ROOT output, several figures and plots are calculated using 
ROOT to confirm the correct execution of the simulation. The benchmark results are 
characterized by: (1) the total number of generated events and detected coincidences; 
(2) their spatial and time distributions; (3) the fractions of random and scattered 
coincidences; and (4) the average acollinearity between the two annihilation gammas. 

The PET benchmark has been run on 12 different system configurations. Two 
operating systems were tested : Linux and Mac OS (versions 10.2.6 and 10.3). The 
Linux distributions were RedHat (versions 9.0, 8.0 and 7.3) and SuSE (version 8.1). 
The source code compilation was performed with either gcc 2.95 or 3.2. The computing 
time for the PET benchmark averaged around 12 and 6 hours for 1.0 GHz and 3 GHz 
processors, respectively. 

PET physical variables characterizing the simulation results are shown in Tabled 
with their mean value and standard deviation obtained from the run of the PET 
benchmark on the 12 system configurations, using a different initial seed on each 
system for the generation of random numbers. The results in Tableland, in particular, 
the relative standard deviations show that the main physical simulation variables are 
stable within less than 1%. It is strongly recommended that the user validates a new 
or updated GATE installation using this table. 
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Table 1. Average values and relative standard deviations (stdev) of the physical 
variables studied with the PET benchmark. 



variable type 



average value relative stdev 



simulated 15-Oxygen life time 
gamma acoUinearity angle 



total decays during the acquisition 
random coincidences 



unscattered coincidences 



scattered coincidences 



3.6815 X 10^ ±0.01% 

20, 568 ±0.58% 

311,778 ±0.39% 

369, 915 ±0.33% 

123.316 S ±0.06% 



0.6063 deg ±0.01% 



5.2. SPECT benchmark 

The SPECT benchmark (Figure |3J) simulates a SPECT acquisition with a moving 
source. The simulated 4-head gamma camera does not correspond to any real system. 
This benchmark involves a cylindrical phantom (5 cm in diameter and 20 cm long) 
filled with water with an inner cylinder (2 cm in diameter, 5 cm long) filled with 
30 kBq of ^^™Tc. The phantom lies on a table (0.6 cm thick, 8 cm wide, and 34 cm 
long). During the simulated acquisition, the table and phantom translate together at 

0. 04 cm/s. Confinement is used to keep the source distribution synchronized with the 
phantom movement. All 4 heads of the gamma camera are identical, consisting of a 
parallel hole lead collimator (hole diameter: 0.3 cm, collimator thickness: 3 cm, and 
septa thickness: 0.6 mm), a 1 cm thick Nal crystal, a 2.5 cm thick back-compartment 
in Perspex, and a 2 cm thick lead shielding. The simulated acquisition consists of 
64 projections (16 projections per head), acquired along a circular orbit with a 7 cm 
radius of rotation and a speed of 0.15 deg/s. Sixteen runs of 37.5 s each are performed 
to simulate the 16 positions of the 4 gamma camera heads. 

The low energy electromagnetic processes package is used to model the physics. 
Rayleigh, photoelectric, and Compton interactions are turned on while the gamma 
conversion interactions are turned off. To speed up the simulation, the X-ray 
production cut is set at 20 keV and secondary electrons are not tracked. Compton 
events occurring in the phantom, collimator, back-compartment, shielding and table 
are recorded. A Gaussian energy blur is applied to all events detected in the crystal, 
using an energy resolution of 10% at 140 keV. The limited spatial resolution of the 
photomultipliers and associated electronics is modeled using a Gaussian blur with a 
standard deviation of 2 mm. Only photons detected with an energy between 20 and 
190 keV are stored. 

The benchmark results are characterized by: (1) the number of simulated events 
and of detected counts; (2) the percentage of primary and scattered events with respect 
to all events detected in the 20-190 keV energy window (here scattered events are 
considered a function of the compartment in which the last scattered event occurred 

1. e. phantom, collimator, table, crystal, or back-compartment); (3) the mean and 
standard deviation of the number of detected counts per projection; (4) the percentage 
of scattered events as a function of the scattering order (1"* for single scatter, 2'"* for 
double scatter, and so on). 

The SPECT benchmark has run on 11 system configurations. Two operating 
systems have been tested: Linux (RedHat versions 9.0, 8.0, 7.3, 7.1, SuSE versions 9.0 
and 8.1, and Fedora Core 1) and Mac OS (version 10.2.8). The mean and standard 
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Figure 5. Illustration of the SPECT benchmark setup. 



Table 2. Average values and relative standard deviations (stdev) of the figures 
of merit used for the SPECT benchmark. 



Global information 






variable type 


average value 


relative stdev 


number of emitted particles 


1.79994 X 10^ 


±0.01% 


detected counts between 20 and 190 keV 


35,919 


±0.3% 


percentage of unscattered photons 


32.9% 


±1.0% 


mean detected counts per projection 


278.4 


±0.9% 


Percentage of photons whose last scattered event occurred in a specific medium 


medium 


average value 


relative stdev 


phantom 


52.3% 


±0.6% 


table 


3.1% 


±1.8% 


collimator 


2.1% 


±2.4% 


crystal 


8.5% 


±1.2% 


back-compartment 


1.2% 


±3.3% 


Percentage of scattered photons as a 


function of the scattering order 


scattering order 


average value 


relative stdev 


order 1 


46.4% 


±0.7% 


order 2 


26.8% 


±0.7% 


order 3 


13.6% 


±1.5% 


order 4 


6.7% 


±0.7% 


order > 4 


6.5% 


±1.6% 



deviation of most figures of merit characterizing the results of the runs are given in 
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Table |2 Similar to the PET benchmark resuhs, this table shows that the results 
produced by GATE are very stable. Table [21 should be used to validate any new 
installation or update of GATE. The time needed to run the benchmark ranged from 
about 3 hours (2.8 GHz Pentium, 2 Gb RAM) to about 11 hours (1 GHz Pentium, 
2.3 Gb RAM). 

6. Validation of GATE 

The validation of Monte Carlo simulated data against real data obtained with 
PET and SPECT cameras is essential to assess the accuracy of GATE and the 
OpenGATE collaboration is largely involved into the validation of GATE. Tables|2|and 
^list commercial systems which have been or are currently being considered for PET 
and SPECT validations. These tables summarize the figures of merit (FOM) used 
for assessing the consistency between simulated and real data, as well as the major 
validation results and associated references. For details regarding these validation 
studies, the reader is highly encouraged to refer to the appropriate references. Overall, 
these studies illustrate the flexibility and reliability of GATE for accurately modelling 
various detector designs. The modelling of the Millennium VG (GE MS) is also in 
progress. 

Table 3. Validation result summary of commercial systems already or currently 
considered for GATE validation in PET. 



PET system 


Studied FOM 


Experiment/GATE 


References 


ECAT EXACT HR+, 


Spatial resolution 




Jan et al 2003 


CPS 


- radial @10cm 


11.4 mm / 11.8 mm 






- tangential @10cm 


11.8 mm / 10.7 mm 






3D sensitivity 


0.75% / 0.80% 






3D count rates 








- prompts @10kBq/ml 


550 kcps / 550 kcps 






- trues OlOkBq/ml 


330 kcps / 300 kcps 






3D scatter fraction 


36% / 35% 




Allegro, Philips 


3D count rate 


@555MBq/ml 


Lamare et al 




- trues 


800 kcps / 950 kcps 


2004 




- scatter 


950 kcps / 900 kcps 






- randoms 


2,000 kcps / 2,400 kcps 






3D scatter fraction 


8% difference 




GE Advance, GE MS 


Energy spectra 


visual assessment 


Schmidtlein et al 




3D scatter fraction 


47.1% / 47.2% 


2004 


MicroPET 4, 


Spatial resolution 




Jan et al 2003b 


Concorde 


- radial (9!2cm 


2.35 mm / 2.25 mm 






- tangential @2cm 


2.45 mm / 2.30 mm 






Sensitivity (350-650 keV) 


1.43% / 2.42% 






Miniature Derenzo 


visual assessment 




MicroPET Focus, 


Spatial resolution 




Jan et al 2004 


Concorde 


- radial @8cm 


3.55 mm / 3.4 mm 






Sensitivity 


3.4% / 3.5% 
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Table 4. Validation result summary of commercial systems already or currently 
considered for GATE validation in SPECT. 



bPEC i system 


Studied h UM 


i!jxperiment/GAi E 


Keierences 


IRIX, Philips 


None reported 


n/a 


Staelens et al 2004b 


AXIS, Philips 


Spatial resolution 


1.30 cm / 1.36 cm 


Staelens et al 2003 




Energy spectra 


visual assessment 






Sensitivity 


231 cps/MBq / 246 cps/MBq 






Scatter profiles 


visual assessment 




Solstice, Philips 


Sensitivity 


good agreement 


Staelens et al 2003 






with theoretical data 


Staelens et al 2004c 








Staelens et al 2004d 


DST Xli, GEMS 


Energy spectra 


excellent agreement 


Assie et al 2004b 




Spatial resolution 








- @10cm in air 


9.5 mm / 9.6 mm 






- @20cm in water 


14.2 mm / 14.4 mm 






Sensitivity 


<4% difference 





GATE has also been shown to be appropriate for simulating various prototype 
imaging devices dedicated to small animal imaging. Table [S] shows the prototypes 
currently simulated using GATE and indicates the features that have been studied 
and validated against experimental measurements. 



Table 5. Prototypes dedicated to small animal imaging modeled using GATE and 
features that have been studied using simulated data and summary validation 
results when available. 



Prototype 


Studied FOM 


Experimcnt/GATE 


Reference 


LSO/LuYAP phoswich PET 


Sensitivity 


n/a 


Rey et al 2003 




NEC curves 


n/a 




High resolution dual head PET 


Spatial resolution 




Chung et al 2003 




- at center 


1.60 mm / 1.55 mm 






- 4 mm off-center 


1.72 mm / 1.7 2mm 






Sensitivity 


0.13% / 0.12% 






Line phantom 


visual assessment 




CsI(Tl) SPECT camera 


Energy spectra 


good agreement 


Lazaro et al 2004 




Spatial resolution 








- @10cm in air 


6.7 mm / 6.8 mm 






Scatter fraction 


0.531% / 0.527% 






Sensitivity 


<2% difl'erence 






Line phantom 


visual assessment 




OPET 


Spatial resolution 


n/a 


Rannou et al 2003 




Sensitivity 


n/a 
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7. Work in progress and future developments 

7.1. Speed-up techniques 

Compared to simpler dedicated codes like SimSET (Harrison et al 1993), the versatility 
of GATE comes at the expense of relatively long computation times. To compensate 
for this, variance reduction tools are currently being developed for GATE. Another 
approach to improve the computing performance of GATE is to distribute the 
simulations on multiple architectures. This is referred to as the gridification of GATE, 
and consists of subdividing simulations on geographically distributed processors in a 
Grid environment by parallelizing the random number generator. GATE simulations 
use a very long period pseudo-random number generator developed from the algorithm 
of F. James (James 1990, Marsaglia and Zaman 1987). This random number generator 
can be subdivided into 900 million different non-overlapping sub-sequences, each 
containing approximately 10^*^ numbers. These sub-sequences provide convenient 
starting points throughout the main sequence. Parallel simulations are produced 
by using the sub-sequences as independent streams by a sequence splitting method 
(Traore and Hill 2001). To demonstrate the potential benefits of gridification, a 
simulation was performed. Table El gives the total computing time in minutes of a 
simulation running on a single 1.5 GHz Pentium IV and the same simulation split 
between 10, 20, 50 and 100 processors on the European DataGrid testbed. In this 
experiment, the testbed consisted of 200 dual-processors with a mix of 750 MHz and 
1 GHz Pentium HI, and 1.4 GHz Pentium IV processors. Gridification using 20 jobs 
decreased the computing time by a factor of 8 when compared to using a single Pentium 
IV 1.5 GHz processor. This example emphasizes the fact that computing time is not 
proportional to the number of jobs running in parallel, due in part to the launch time 
of the jobs, the modelling of the geometry of the simulation (which is independent 
of random number generation), and the time spent in queuing jobs. This study also 
proved that the results obtained via the Grid were equivalent to those generated on a 
single machine. In the future, each Unix platform on the DataGrid will be installed 
with GATE, allowing the use of more than 500 processors for a simulation. The 
development of a convenient tool to split, launch, and retrieve GATE simulations on 
a Grid environment using a web interface is currently under development. 

Table 6. Comparison of computing times between local and parallel jobs. 
Total computing time in minutes 



1 1.5 GHz Pentium IV 159 

10 jobs 31 

20 jobs 21 

50 jobs 31 

100 jobs 38 



7.2. Extension of GATE to other domains 

Besides its application to classical emission tomography, GATE is potentially appro- 
priate to simulate in-line tomography in hadrontherapy. Indeed, GATE possesses 
the ability to model the distribution of Z?"*" emitters along the beam path. In prin- 
ciple, this distribution can be obtained from knowledge of the nuclear cross-sections 
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T GATE 




Angle 6 (degrees) 

Figure 6. Anisotropy functions at r = 2 cm for a ^^'^Ir brachytherapy source. 
The anisotropy functions F{r,t) have been normahzed to one at 90 deg. 



of the reactions between heavy ions and target nuclei. This distribution can then be 
accurately reproduced in GATE using a voxelized source, composed of the different 
emitters PO, "C and ^°C). 

Currently, Monte Carlo simulations are believed to be the most accurate method 
for dose calculation in radiotherapy and brachytherapy. A comparison of the 
computational codes in radiation dosimetry is under way (Berger et al n.d.). In this 
context, the anisotropy function F{r, 9) (Nath et al 1995) of a i^^Jr brachytherapy 
source in water has been studied. As illustrated in Figure El anisotropy functions in 
water have been calculated using Geant3 (Brun et al 1987), MCNP4C (Briesmeister 
1993), MCPT (WiUiamson 1998), and GATE at a radial distance of 2 cm. GATE was 
in agreement with all other codes. The relative deviation concerning anisotropy 
functions between MCPT and GATE simulations is less than 3%. This indicates that 
GATE has potential to perform dose calculations in brachytherapy. The challenge 
for GATE is to provide detailed descriptions of local doses delivered by radiotherapy 
or brachytherapy treatment. It is envisioned that this can be done with voxelized 
phantoms in future GATE simulations. 

8. Conclusion 

Based on Geant4, GATE is a versatile and adaptable platform for simulating PET 
and SPECT experiments. GATE is appropriate for simulating conventional scanners 
and novel detection devices, and does not require any knowledge of CH — h. The 
OpenGATE collaboration, representing a large number of research groups from around 
the world, has publicly released this simulation toolkit after two years of software 
development and validation. The source code is available for download, which will 
enable users to make modifications to suit their particular needs. Documentation is 
available and user support is provided through a very active mailing list. 
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The future of GATE is closely related to the future of Monte Carlo simulations in 
nuclear medicine. Because Monte Carlo simulations are playing an increasing role in 
the optimization of detector design and in the assessment of acquisition and processing 
protocols, we hope that GATE will answer many of these needs. In addition, the efforts 
of the OpenGATE collaboration should ensure that GATE will continue to evolve 
to become a comprehensive simulation tool at the service of the nuclear medicine 
community. 
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